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the  differences  in  source  function. 


High  frequency  phase  ratios  offer  some  hope  of  discrimination,  but  at  the  cost  of 
increased  scatter  and  noise.  Also,  the  California  events  show  a  distinct  trend  of  increasing 
Pg/Lg  with  distance.  This  trend  is  accentuated  at  high  frequencies,  thus  requiring  a  distance 
correction  for  effective  discrimination. 
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INTRODUCTION 


Some  of  the  more  promising  means  of  discriminating  between  explosions  and  earth¬ 
quakes  involve  the  ratio  of  P  waves  to  Lg  or  S  waves.  For  example,  Dysart  and  Pulli  (19X7) 
used  Pn/Sn  and  Pn/Lg  ratios  to  discriminate  between  regional  earthquakes  and  mine  blasts 
recorded  at  the  NORRSS  array.  With  only  a  few  exceptions,  explosions  had  higher  Pn/Lg 
ratios  than  earthquakes.  The  explosions  also  tended  to  have  higher  Pn/Sn  ratios  than  the 
earthquakes.  On  the  other  hand,  similar  studies  by  Hutchcnson  and  Kraft  (1986)  and  Bennett 
and  Murphy  (1986)  reported  somewhat  mediocre  discrimination  of  explosions  and  earthquakes 
in  the  western  U.  S.  by  P/Lg  ratios.  Taylor  et  at.  (1989)  confirmed  the  latter  results  using  the 
Lawrence  Livermore  National  Laboratory  (LLNL)  broadband  network.  The  reasons  for  the 
uneven  performance  of  the  P/Lg  discriminant  are  unclear,  but  they  are  vital  for  the  evaluation 
and  application  of  the  discriminant.  A  comprehensive  review  of  early  studies  on  P/Lg  is 
given  in  Pomeroy  et  al.  (1982). 

Several  factors  may  contribute  to  this  uneven  performance.  The  most  obvious  is  the  tec¬ 
tonic  region:  NORESS  is  located  in  the  Precambrian  Baltic  shield  area,  and  most  of  the 
sources  in  the  Dysart  and  Pulli  (1987)  study  occurred  in  the  shield  or  in  Paleozoic  orogenic 
terranes.  On  the  other  other  hand,  the  western  U.  S.  studies  involve  source-station  geometries 
located  largely  in  the  Tertiary  basin-and-range  and  orogenic  terranes.  The  large  differences  in 
crustal  thickness,  velocity  structure  and  attenuation  may  be  responsible  for  the  differences 
between  the  studies.  Another  complication  is  the  frequency:  Hutchenson  and  Kraft  (1986) 
measured  P/Lg  ratios  for  frequencies  around  1  Hz,  whereas  the  NORESS  results  used  higher 
frequencies  (3-12  Hz).  Furthermore,  a  later  study  by  Dysart  and  Pulli  (1988)  of  NORESS 
data  using  intermediate  frequencies  was  unsuccessful  in  discriminating  earthquakes  and 
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explosions.  On  the  other  hand,  Taylor  el  al.  (1989)  found  no  significant  difference  between 
~1  Hz  and  6-8  Hz  data  for  the  LLNL  broadband  network  in  the  western  U.  S. 

One  of  the  difficulties  in  determining  why  the  discriminant  fails  is  that  the  reasons  for  its 
success  are  still  problematic.  One  obvious  possibility  is  simply  that  explosions  generate  less 
shear  energy  than  earthquakes  for  a  given  source  size:  the  main  sources  of  shear  energy  from 
explosions  are  P— >S  conversion  at  the  free  surface  and  scattering.  Another  possibility,  pro¬ 
posed  by  Frankel  (1989),  is  that  the  P/Lg  ratio  is  a  depth  discriminant  due  to  very'  high 
attenuation  near  the  surface,  where  the  explosions  are  detonated.  Another  possible  explana¬ 
tion  for  the  discriminant’s  apparent  success  is  that  the  explosions  in  these  studies  generally 
come  from  a  fairly  limited  area,  which  is  often  different  from  the  source  areas  of  the  earth¬ 
quakes,  suggesting  that  the  difference  in  P/Lg  ratio  between  the  two  populations  is  due  simply 
to  path  effects  or  near-source  geology. 

In  this  study,  we  have  attempted  to  identify  the  reasons  for  the  success/failure  of  the 
P/Lg  discriminant  in  the  western  U.  S.  using  seismograms  recorded  at  the  DWWSSN  station 
JAS  (Jamestown,  California).  We  have  selected  earthquakes  whose  source-station  paths  are  as 
similar  as  possible  to  those  for  the  explosions,  all  of  which  were  detonated  at  NTS  (Figure  1). 
The  events  used  are  listed  in  Tables  1  and  2.  We  have  also  studied  a  wide  range  of  frequen¬ 
cies  in  order  to  determine  whether  a  particular  frequency  band  provides  significantly  better  or 
worse  performance  than  others.  Also,  three  pairs  of  explosions  and  their  resultant  collapses 
have  been  studied  to  eliminate  the  effect  of  propagation  path  on  P/Lg  ratios. 

The  primary  results  from  this  can  be  summarized  as  follows: 

(1)  The  Pg/Lg  ratio  for  explosions  varies  strongly  as  a  function  of  source  depth  and 
size,  so  that  smaller  sources  appear  more  nearly  similar  to  earthquakes  than  do  the 
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larger  sources. 


(2)  The  explosion  Pg/Lg  ratios  are  actually  similar  to  earthquake  Pg/Lg  ratios  along  the 
same  travel  pith,  whereas  both  are  quite  different  from  earthquake  Pg/Lg  ratios 
along  a  significantly  different  path.  This  suggests  that  path  effects  are  at  least  as 
important  as  source  type  in  determining  Pg/Lg  ratios  in  the  western  U.  S. 

(3)  In  the  highest  frequency  range  studied,  4-8  Hz,  the  path  effects  seem  to  be  dimin¬ 
ished,  and  separation  between  explosions  and  earthquakes  is  improved  over  that  in 
the  lower  frequency  bands. 
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Table  1. 

Explosions  used  in  this  study. 

Date 

Time 

Latitude 

Longitude 

Depth 

Magnitude 

Name 

14  NOV  1980 

16:50:00 

37.11 

-116.02 

0.320 

4.5 

DAUPHIN 

17  DEC  1980 

15:10:00 

37.30 

-116.30 

0.573 

5.0 

SERPA 

29  MAY  1981 

16:00:00 

37.10 

-116.00 

0.320 

4.5 

ALIGOTE 

10  JUL  1981 

14:00:00 

37.13 

-116.03 

0.340 

4.2 

NIZA 

16  JUL  1981 

15:00:00 

37.09 

-116.02 

0.204 

3.3 

PINEAU 

27  AUG  1981 

14:31:00 

37.16 

-116.07 

0.294 

4.3 

ISLAY 

24  SEP  1981 

15:00:00 

37.01 

-116.02 

0.213 

3.5 

CERNADA 

1  OCT  1981 

19:00:00 

37.10 

-116.00 

0.472 

5.2 

PALIZA 

11  NOV  1981 

20:00:00 

37.08 

-116.07 

0.445 

5.0 

TILCI 

12  NOV  1981 

15:00:00 

37.10 

-116.00 

0.518 

5.3 

ROUSANNE 

3  DEC  1981 

15:00:00 

37.15 

-116.07 

0.494 

4.9 

AKAVI 

16  DEC  1981 

21:05:00 

37.12 

-116.12 

0.335 

4.4 

CABOC 

28  JAN  1982 

16:00:00 

37.10 

-116.10 

0.640 

5.8 

JORNADA 

29  JUL  1982 

20:05:00 

37.10 

-116.08 

0.400 

4.6 

MONTEREY 

5  AUG  1982 

14:00:00 

37.10 

-116.00 

0.640 

5.4 

ATRISCO 

2  SEP  1982 

14:00:00 

37.02 

-116.02 

0.229 

3.3 

CERRO 

23  SEP  1982 

16:00:00 

37.20 

-116.20 

0.408 

4.8 

HURONLANDING/ 

DIAMONDACE 

23  SEP  1982 

17:00:00 

37.20 

-116.10 

0.451 

4.8 

FRISCO 

29  SEP  1982 

13:30:00 

37.09 

-116.05 

0.564 

4.1 

BORREGO 

12  NOV  1982 

19:17:00 

37.02 

-116.03 

0.366 

4.3 

SEYVAL 

1 1  FEB  1983 

16:00:00 

37.10 

-116.00 

0.304 

4.1 

COALORA 

17  FEB  1983 

17:00:00 

37.16 

-116.06 

0.343 

4.0 

CHEEDAM 

26  MAY  1983 

15:00:00 

37.10 

-116.00 

0.384 

4.4 

FAHADA 

9  JUN  1983 

17:10:00 

37.16 

-116.09 

0.320 

4.6 

DANABLU 

3  AUG  1983 

13:33:00 

37.12 

-116.09 

0.326 

4.3 

LABAN 

11  AUG  1983 

14:00:00 

37.00 

-116.00 

0.320 

4.2 

SAB ADO 

22  SEP  1983 

15:00:00 

37.11 

-116.05 

0.533 

4.1 

TECHADO 

31  JAN  1984 

15:30:00 

37.11 

-116.12 

0.388 

4.4 

GORBEA 

20  JUN  1984 

15:15:00 

37.00 

-116.00 

0.381 

4.5 

DUORO 

2  AUG  1984 

15:00:00 

37.00 

-116.00 

0.335 

4.4 

CORREO 

30  AUG  1984 

14:45:00 

37.09 

-116.00 

0.366 

4.3 

DOLCETTO 

13  SEP  1984 

14:00:00 

37.10 

-116.10 

0.483 

5.0 

BRETON 

23  MAR  1985 

18:30:00 

37.18 

-116.09 

0.515 

5.1 

COTTAGE 

6  APR  1985 

23:15:00 

37.20 

-116.20 

0.000 

4.8 

MISTYRAIN 

12  JUN  1985 

17:30:00 

37.10 

-116.10 

0.293 

4.2 

VILLE 

26  JUN  1985 

18:03:00 

37.12 

-116.12 

0.381 

4.0 

MARIBO 

17  AUC.  1985 

16:25:00 

37.00 

-116.04 

0.332 

4.2 

CHAMITA 

16  OCT  1985 

21:35:00 

37.11 

-116.12 

0.415 

4.6 

ROQUEFORT 

22  MAR  1986 

16:15:00 

37.10 

-116.10 

0.610 

5.0 

GLENCOE 

10  APR  1986 

14:08:30 

37.20 

-1 16.20 

0.400 

4.9 

MIGHTYOAK 

21  MAY  1986 

13:59:00 

37.13 

-116.06 

0.500 

3.9 

PANAMINT 

5  JUN  1986 

15:04:00 

37.10 

-1 16.02 

0.500 

5.3 

TAJO 

24  JUL  1986 

15:05:00 

37.10 

-116.10 

0.400 

4.4 

CORNUCOPIA 

5 


Table  2.  Earthquakes  used  in  this  study. 

Date 

Time 

Latitude 

Longitude 

Depth 

Magnitude 

2  OCT  1980 

01:48:13 

37.30 

-117.00 

5 

3.0 

18  OCT  1980 

03:54:30 

34.40 

-116.70 

8 

3.4 

19  OCT  1980 

22:59:52 

34.40 

-116.70 

9 

3.1 

13  DEC  1980 

04:10:16 

35.90 

-118.00 

5 

3.1 

13  DEC  1980 

14:15:47 

34.50 

-116.30 

5 

3.4 

29  DEC  1980 

07:12:53 

37.50 

-113.10 

7 

3.1 

2  MAY  1981 

02:11:05 

35.90 

-117.80 

5 

3.2 

5  MAY  1981 

13:59:05 

36.40 

-118.10 

1 

3.8 

5  MAY  1981 

14:34:52 

36.40 

-118.00 

1 

3.7 

18  MAY  1981 

21:32:48 

37.20 

-117.90 

6 

3.1 

23  MAY  1981 

18:50-21 

36.10 

-117.90 

1 

3.1 

11  JUN  1981 

18:00:43 

38.30 

-115.90 

5 

3.6 

5  JUL  1981 

10:30:48 

35.80 

-117.70 

5 

3.1 

5  JUL  1981 

12:31:00 

35.80 

-117.70 

5 

3.0 

5  JUL  1981 

13:33:50 

35.80 

-117.70 

5 

3.0 

5  JUL  1981 

17:45:24 

35.80 

-117.70 

7 

3.0 

25  JUL  1981 

20:05:54 

36.10 

-117.80 

6 

3.0 

23  AUG  1981 

04:32:29 

38.00 

-117.60 

6 

3.0 

12  SEP  1981 

21:23:07 

34.20 

-117.30 

7 

3.6 

16  SEP  1981 

07:57:42 

34.20 

-117.30 

4 

3.0 

13  OCT  1981 

14:47:53 

37.10 

-116.90 

6 

3.4 

13  OCT  1981 

19:51:12 

37.00 

-117.00 

5 

3.1 

15  OCT  1981 

04:21:10 

37.10 

-116.90 

1 

3.5 

11  NOV  1981 

20:25:10 

35.80 

-117.70 

7 

3.2 

19  NOV  1981 

21:40:53 

37.10 

-116.90 

6 

3.1 

24  NOV  1981 

09:20:17 

34.80 

-116.70 

6 

3.5 

16  DEC  1981 

01:34:52 

36.20 

-117.90 

3 

3.5 

24  JAN  1982 

15:44:07 

37.50 

-117.80 

5 

4.1 

28  JAN  1982 

22:51:02 

38.50 

-118.10 

5 

4.5 

31  JUL  1982 

00:57:58 

35.80 

-117.70 

5 

3.2 

14  AUG  1982 

02:37:59 

34.20 

-117.30 

5 

3.1 

2  SEP  1982 

22:51:06 

38.10 

-117.90 

5 

3.2 

25  SEP  1982 

23:29:46 

36.40 

-117.80 

5 

3.5 

28  SEP  1982 

10:43:51 

35.80 

-117.80 

5 

3.5 

28  SEP  1982 

17:35:05 

35.80 

-117.80 

9 

3.1 

29  SEP  1982 

18:19:16 

35.80 

-117.80 

6 

3.9 

29  SEP  1982 

18:21:01 

35.80 

-117.80 

8 

4.2 

29  SEP  1982 

19:37:14 

35.80 

-117.80 

9 

3.8 

30  SEP  1982 

22:38:10 

35.80 

-117.80 

8 

4.1 

1  OCT  1982 

01:33:35 

37.90 

-118.20 

5 

3.9 

1  OCT  1982 

06:24:58 

35.80 

-117.80 

4 

3.1 

1  OCT  1982 

09:21:16 

35.80 

-117.80 

4 

3.2 

1  OCT  1982 

12:19:32 

35.70 

-117.80 

5 

3.4 

1  OCT  1982 

17:14:42 

35.80 

-117.80 

8 

3.8 

6 


Table  2.  Earthquakes  used  in  this  study. 

Dale 

Time 

Latitude 

Longitude 

Depth 

Magnitude 

1  OCT  1982 

20:01:26 

35.80 

-117.80 

8 

3.1 

1  OCT  1982 

20:45:55 

35.80 

-117.70 

6 

3.7 

1  OCT  1982 

20:46:03 

35.80 

-117.70 

6 

3.8 

1  OCT  1982 

21:14:13 

35.70 

-117.70 

5 

3.1 

1  OCT  1982 

22:10:21 

35.70 

-117.80 

7 

4.3 

2  OCT  1982 

09:33:05 

37.90 

-118.30 

5 

3.7 

2  OCT  1982 

14:01:56 

35.80 

-117.70 

6 

3.5 

2  OCT  1982 

16:01:21 

35.80 

-117.80 

8 

3.5 

3  OCT  1982 

09:47:42 

35.80 

-117.80 

5 

3.0 

4  OCT  1982 

02:23:59 

37.90 

-118.10 

5 

j>.9 

4  OCT  1982 

18:43:28 

35.80 

-117.80 

8 

3.9 

6  OCT  1982 

11:38:40 

35.80 

-117.60 

4 

3.0 

7  OCT  1982 

15:50:05 

35.70 

-117.70 

5 

3.2 

7  OCT  1982 

17:54:36 

35.70 

-117.70 

10 

3.9 

12  OCT  1982 

08:22:46 

35.80 

-117.70 

8 

3.7 

12  OCT  1982 

17:32:31 

35.80 

-117.70 

9 

3.1 

21  OCT  1982 

19:23:35 

35.70 

-117.70 

6 

3.3 

29  OCT  1982 

09:20:20 

35.70 

-117.70 

6 

3.2 

30  OCT  1982 

14:43:13 

35.70 

-117.80 

7 

3.1 

10  NOV  1982 

11:21:25 

34.00 

-116.70 

8 

3.9 

12  NOV  1982 

05:26:47 

38.00 

-118.10 

5 

3.9 

1  JAN  1983 

04:12:21 

35.80 

-117.70 

6 

3.0 

2  JAN  1983 

16:32:19 

36.50 

-116.60 

5 

3.2 

3  JAN  1983 

17:39:43 

36.50 

-116.60 

5 

3.1 

5  JAN  1983 

09:17:16 

37.60 

-117.80 

6 

3.2 

5  JAN  1983 

16:06:34 

35.70 

-117.80 

8 

3.0 

6  JAN  1983 

20:30:51 

35.80 

-117.70 

7 

3.0 

7  JAN  1983 

13:43:43 

35.70 

-117.80 

7 

3.9 

7  JAN  1983 

14:03:22 

35.80 

-117.70 

6 

3.0 

7  JAN  1983 

14:19:55 

35.70 

-117.80 

9 

3.0 

7  JAN  1983 

18:15:32 

35.80 

-117.80 

6 

3.1 

7  JAN  1983 

19:22:32 

35.80 

-117.80 

9 

3.0 

8  JAN  1983 

07:19:30 

34.10 

-117.40 

5 

4.1 

8  JAN  1983 

07:30:58 

35.70 

-117.80 

9 

3.3 

8  JAN  1983 

13:26:30 

35.80 

-117.70 

6 

3.1 

13  JAN  1983 

15:00:59 

35.70 

-117.80 

9 

3.0 

17  JAN  1983 

01:09:00 

35.70 

-118.60 

7 

3.7 

1  FEB  1983 

22:09:07 

35.80 

-117.70 

11 

3.1 

23  FEB  1983 

11:10:20 

36.00 

-114.70 

5 

3.9 

1  MAR  1983 

01:49:16 

35.80 

-117.70 

6 

3.3 

2  MAR  1983 

14:23:55 

35.80 

-117.70 

10 

3.2 

4  MAR  1983 

05:51:20 

35.80 

-117.70 

10 

3.0 

1 1  MAR  1983 

05:26:29 

35.80 

-117.70 

10 

3.1 

21  MAY  1983 

11:02:31 

36.10 

-117.10 

5 

3.4 

7 


Table  2.  Earthquakes  used  in  this  study. 

Dale 

Time 

Latitude 

Longitude 

Depth 

Magnitude 

4  JUN  1983 

11:37:40 

37.40 

-115.20 

6 

3.6 

15  JUL  1983 

22:03:10 

38.50 

-118.30 

13 

4.6 

18  JUL  1983 

04:36:47 

34.20 

-117.10 

3 

3.2 

12  SEP  1983 

12:08:02 

34.00 

-117.30 

15 

3.6 

12  SEP  1983 

21:33:35 

35.20 

-117.00 

6 

3.2 

23  SEP  1983 

01:35:26 

35.60 

-118.60 

2 

3.2 

19  OCT  1983 

14:00:37 

35.90 

-118.30 

5 

4.2 

19  OCT  1983 

19:48:24 

35.90 

-118.30 

5 

3.8 

30  OCT  1983 

20:02:59 

35.90 

-118.30 

5 

3.9 

6  NOV  1983 

01:34:31 

35.90 

-118.30 

1 

3.1 

6  NOV  1983 

02:08:1' 

35.90 

-118.30 

0 

3.1 

9  DEC  1983 

08:58:41 

38.60 

-112.60 

7 

3.6 

18  DEC  1983 

17:56:48 

35.80 

-118.00 

4 

3.3 

29  DEC  1983 

07:23:29 

35.90 

-118.30 

1 

3.3 

29  DEC  1983 

19:46:16 

34.20 

-117.40 

8 

3.6 

29  DEC  1983 

19:49:13 

34.20 

-117.40 

6 

2.9 

29  DEC  1983 

19:50:02 

34.20 

-117.30 

8 

2.7 

12  JAN  1984 

15:47:16 

35.90 

-118.30 

1 

3.1 

16  JAN  1984 

04:06:12 

35.90 

-118.30 

0 

3.4 

20  JAN  1984 

08:38:49 

36.20 

-118.30 

6 

3.1 

20  JAN  1984 

12:07:50 

36.00 

-118.00 

9 

3.3 

11  JUN  1984 

22:21:10 

34.40 

-116.60 

2 

4.0 

17  JUN  1984 

20:23:24 

38.20 

-117.90 

6 

3.8 

19  AUG  1984 

02:22:34 

38.40 

-118.10 

5 

4.2 

7  SEP  1984 

14:51:51 

36.00 

-117.30 

0 

3.7 

2  JAN  1985 

05:24:58 

34.00 

-116.50 

9 

3.8 

8  JAN  1985 

06:59:38 

36.90 

-118.00 

5 

3.6 

1  APR  1985 

06:13:33 

36.00 

-117.40 

0 

3.2 

13  MAY  1985 

21:24:00 

35.80 

-117.70 

6 

3.4 

10  JUN  1985 

00:58:01 

34.20 

-116.80 

11 

3.3 

18  JUN  1985 

01:23:40 

35.20 

-117.30 

8 

3.8 

16  JUL  1985 

17:57:50 

34.50 

-116.80 

0 

3.9 

6  AUG  1985 

03:45:36 

35.40 

-117.70 

12 

3.1 

14  AUG  1985 

06:12:56 

35.00 

-116.90 

8 

3.6 

16  AUG  1985 

01:51:21 

36.20 

-117.90 

5 

4.3 

19  AUG  1985 

23:10:26 

34.90 

-116.90 

6 

3.2 

29  AUG  1985 

07:59:08 

34.30 

-116.80 

5 

4.0 

2  SEP  1985 

19:45:08 

37.20 

-118.00 

5 

3.4 

10  DEC  1985 

06:10:25 

37.30 

-115.00 

5 

3.7 

16  DEC  1985 

09:09:52 

36.30 

-118.00 

6 

3.2 

2  MAR  1986 

23:03:01 

35.90 

-118.30 

6 

3.2 

6  MAR  1986 

20:16:52 

37.20 

-117.30 

5 

3.7 

21  APR  1986 

06:35:59 

35.80 

-117.80 

4 

3.3 

23  MAY  1986  11:41:55  35.80  -118.00 _ 10 _ ^9 


DATA  AND  RESULTS 


Methodology 

The  first  station  studied  was  JAS  (Jamestown,  California),  due  to  its  long  digital  history 
and  large  number  of  recorded  events  from  the  Nevada  Test  Site  and  the  California-Nevada 
border  region.  The  station  was  moved  in  1984  and  renamed  JAS1,  but  since  this  later  site  is 
less  than  3  km  away,  both  data  sets  are  considered  together.  The  events  were  selected  in 
order  to  minimize  the  differences  in  epicentral  distance  and  azimuth  in  an  effort  to  limit  the 
effect  of  gross  path  differences  (Figure  1).  Since  well  recorded  explosions  in  this  area  were 
basically  limited  to  the  Nevada  Test  site,  all  such  explosions  were  used.  Earthquakes  were 
chosen  to  be  as  close  to  NTS  as  possible. 

The  seismograms  were  processed  by  picking  the  onset  of  the  Pg  and  Lg  phases  and 
measuring  the  Pg/Lg  ratio  in  the  frequency  bands  0.5-1.0  Hz,  1-2  Hz,  2-4  Hz,  3-6  Hz,  and  4-8 
Hz.  An  example  of  both  earthquake  and  explosion  seismograms  is  shown  in  Figure  2, 
bandpassed  in  several  of  these  frequency  bands.  The  explosion  actually  generates  substantial 
low-frequency  Lg  energy,  resulting  in  a  Pg/Lg  ratio  that  is  less  than  that  of  the  earthquake  in 
the  0.5- 1.0  Hz  band.  As  the  frequency  increases,  however,  the  explosion  Lg  energy  drops 
precipitously,  resulting  in  a  Pg/Lg  ratio  that  is  higher  than  that  of  the  earthquake  in  the 
highest  frequency  band. 

One  spectral-domain  and  three  time-domain  measurement  methods  were  tried,  in  an 
effort  to  find  the  most  stable  and/or  effective  method  of  computing  phase  ratios.  In  the 
spectral-domain  method,  12.8  second  windows  were  extracted,  starting  1.3  seconds  before  the 
arrival  onset,  and  cosine  tapers  were  applied  to  the  first  10%  and  last  10%  of  the  window. 
Noise  windows  before  the  Pg  onset  and  immediately  prior  to  the  Lg  onset  were  also  extracted. 
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The  Lg  "noise"  window  actually  consists  largely  of  Pg  coda.  Spectral  amplitudes  were 
obtained  Oy  applying  a  Hartley  transform  (Bracewell,  1986)  to  the  windows,  and  instrument 
corrections  were  applied.  A  geometric  mean  over  the  frequency  band  of  interest  was  then 
computed. 

For  the  time-domain  measurements,  the  seismograms  were  filtered  with  a  zero-phase 
six-pole  Butterworth  bandpass  filter  over  the  frequency  band  of  interest.  Three  types  of  time- 
domain  measurements  were  then  obtained  from  the  10  second  windows  starting  at  the  arrival 
onset: 

(1)  the  root-mean-square  (RMS)  of  the  window 

(2)  the  maximum  zero-to-peak  amplitude  within  the  window,  and 

(3)  the  range  of  the  maximum  peak  to  maximum  trough  within  the  window. 

The  latter  two  time-domain  methods  proved  to  be  somewhat  unstable  at  high  frequencies: 
even  small  glitches  in  the  seismograms,  difficult  to  remove  with  conventional  despiking  algo¬ 
rithms,  produce  anomalously  large  amplitudes  in  the  higher  passbands.  These  glitches  also 
affect  the  spectral  and  the  time-domain  RMS  methods,  but  the  effects  are  mitigated  by  the 
averaging  nature  of  these  methods.  In  the  large,  the  results  were  similar  to  the  spectral 
domain  method,  although  the  latter  two  time-domain  methods  seemed  to  show  somewhat 
more  scatter.  (The  RMS  method  is  analogous  to  the  spectral  method  without  instrument 
correction.)  Throughout  the  rest  of  the  paper,  the  results  will  be  those  obtained  with  the 
spectral-domain  method. 

The  noise  windows  are  important  for  assessing  the  quality  of  a  PgILg  ratio.  In  cases 
where  the  Pg  signal/noise  amplitude  ratio  is  less  than  1.4  (corresponding  to  a  signal/noise 
power  ratio  of  2),  while  the  Lg  signal/noise  amplitude  ratio  is  greater  than  1.4,  the  PgILg  ratio 
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is  actually  an  upper  bound.  In  cases  where  the  Pg  signal/noise  is  greater  than  1.4,  while  the 
Lg  signal/noise  is  less  than  1.4,  the  Pg/Lg  ratio  is  a  lower  bound.  (In  cases  where  both 
signal/noise  ratios  were  less  than  1.4,  the  datum  was  discarded.)  A  lower-bound  value  can  be 
used  in  discrimination  if  it  falls  above  the  discrimination  threshold,  and  an  upper-bound  value 
can  be  useful  if  it  falls  below  the  threshold. 


2  Oct  1980  Earthquake 


16  Jul  1981  Explosion  (PINEAU) 


M  =  3.0  A  =  2.8 


M  =  3.3  A  =  3.6 


0  40  80 

Time  (sec) 


0  40  80 

Time  (sec) 


Figure  2.  Example  seismograms  from  an  earthquake  and  explosion  of  similar  size  and  dis¬ 
tance  from  JAS.  The  top  trace  is  unfiltered;  the  lower  traces  are  filtered  using  a  six-pole 
zero-phase  Butterworth  bandpass  filter. 
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Variation  of  Pg/Lg  with  Distance 


The  Pg/Lg  ratios  for  all  events  in  four  frequency  bands  are  shown  plotted  against 
source-station  distance  in  Figures  3  and  4.  In  all  of  the  frequency  bands,  the  average  Pg/Lg 
ratio  for  explosions  is  higher  than  that  of  the  earthquakes.  However,  the  explosion  and  earth¬ 
quake  populations  exhibit  an  extensive  overlap  in  the  lowest  frequency  band.  As  the  fre¬ 
quency  increases,  the  overlap  diminishes  somewhat,  so  that  the  two  populations  do  show  a 
fair  degree  of  separation  in  the  highest  frequency  band.  Despite  the  still  incomplete  separa¬ 
tion,  the  Pg/Lg  ratio  in  the  4-8  Hz  band  is  nevertheless  useful  in  that  only  5%  of  the  explo¬ 
sions  have  Pg/Lg  ratios  below  0.8,  while  only  2%  of  the  earthquakes  have  Pg/Lg  ratios  above 
1.5.  About  14  of  the  events  fall  into  the  intermediate  range. 

One  possible  explanation  of  the  overlap  is  a  simple  obfuscation  of  explosion  Lg  by  the 
Pg  coda.  In  general,  however,  both  the  low  explosion  Pg/Lg  ratios  and  the  high  earthquake 
Pg/Lg  ratios  tend  to  be  uncensored  values  (i.e.,  both  Pg  and  Lg  have  a  signal/noise  ratio 
greater  than  1.4). 

Due  to  the  different  propagation  characteristics  of  Pg  and  Lg,  including  both  attenuation 
and  geometric  spreading,  it  is  reasonable  to  expect  some  trends  with  distance.  The  scatter  of 
the  data  is  generally  too  great  to  discern  well  defined  trends.  However,  as  the  frequency 
increases,  the  Pg/Lg  ratio  does  appear  to  show  a  general  increase  with  distance.  As  will  be 
shown  later,  this  trend  is  accentuated  by  separating  the  ratios  according  to  the  source-station 
path. 
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0. 5-1.0  Hz 


- 

o 

Earthquakes 

* 

E  xplosions 

- 

V 

upper  bound 

upper  bound 

* 

A 

lower  bound 

* 

lower  bound 

* 


J _ I _ I _ I _ L 

2  3  4  5  6 

Distance  (deg) 


1-2  Hz 


Distance  (deg) 


Figure  3.  Pg/Lg  ratios  plotted  vs.  epicentral  distance  for  the  0.5-1. 0  Hz  and  1-2  Hz  bands. 
Explosions  are  shown  as  crosses  (uncensored  values)  and  arrows  (upper  and  lower  bound 
values).  Earthquakes  are  shown  as  diamonds  (uncensored  values)  and  triangles  (upper  and 
lower  bound  values). 
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Figure  4.  Same  as  Figure  3  for  2-4  Hz  and  4-8  Hz  bands. 
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Variation  of  Pg/Lg  with  Magnitude  and  Depth 


The  Pg/Lg  ratios  are  plotted  against  the  source  magnitude  in  Figures  5  and  6.  The  mag¬ 
nitudes  are  NE1C  body-wave  magnitudes,  if  available.  If  they  are  not  available,  local  magni¬ 
tudes  are  shown.  For  all  frequency  bands,  the  explosion  Pg/Lg  ratios  increase  with  magni¬ 
tude,  whereas  the  earthquake  Pg/Lg  ratios  show  no  apparent  trend  in  any  frequency  band. 
The  result  is  that  small  explosions  tend  to  resemble  earthquakes  insofar  as  their  Pg/Lg  ratios 
are  concerned.  In  fact,  part  of  the  separation  in  explosion  and  earthquake  populations  in  Fig¬ 
ures  3  and  4  is  due  to  a  magnitude  bias  between  the  two  groups. 

The  trend  of  increasing  explosion  Pg/Lg  ratios  with  magnitude  was  also  noted  by  Taylor 
et  al.  (1989).  As  one  might  expect  from  the  strong  correlation  between  explosion  size  and 
depth,  Pg/Lg  for  explosions  also  increases  with  increasing  depth  (Figures  7  and  8)  in  three  out 
of  the  four  frequency  bands.  This  was  also  noted  by  Taylor  et  al.  for  -1  Hz  data  for  the 
LLNL  network,  who  observed  that  the  trend  seemed  to  disappear  at  depths  below  1  km. 
There  does  not  appear  to  be  a  similar  relationship  between  Pg/Lg  and  depth  for  earthquakes 
(Figures  9  and  10).  However,  the  depths,  which  are  obtained  from  NEIC,  are  probably  not 
reliable  enough  to  conclude  that  there  is  no  relationship  between  Pg/Lg  and  depth  for  earth¬ 
quakes. 

One  possible  explanation  for  the  increase  of  Pg/Lg  with  depth  might  be  pP  interference 
with  Pg.  However,  the  shallower  explosions  in  particular  should  be  unaffected  by  pP  interfer¬ 
ence  in  the  0.5- 1.0  Hz  band.  Another  possibility  is  differential  contributions  from  S*  to  the 
Pg  and  Lg  wavetrains,  but  this  effect  should  be  minimal  at  high  frequencies.  Thus,  the 
former  explanation  may  explain  the  trend  in  the  high  frequency  bands,  but  not  the  lowest, 
while  the  latter  may  explain  the  lowest  band  but  not  the  higher.  Perhaps  the  simplest  expla- 
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nation  is  that  the  shallower  explosions  are  detonated  in  lower  velocity  media.  Lower  veloci¬ 
ties  constrict  the  cone  of  rays  contributing  to  Pg  and  thus  the  proportion  of  energy  radiated  as 
Pg.  Furthermore,  Lilwall  (1988)  and  Frankel  (1989)  have  shown  that  the  lower  velocities 
near  the  surface  also  tend  to  increase  Lg  generation,  an  effect  which  increases  with  increasing 
frequency.  In  any  case,  the  source  depth  (or  magnitude)  appears  to  be  a  crucial  factor  in 
determining  the  Pg/Lg  ratio. 
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Figure  6.  Same  as  Figure  5  for  the  2-4  Hz  and  4-8  Hz  bands. 
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Figure  7.  Explosion  Pg/Lg  ratios  for  the  0.5-1. 0  Hz  and  1-2  Hz  band  plotted  vs.  burial  depth 
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Figure  9.  Pg/Lg  ratios  for  all  events  plotted  vs.  source  depth  for  the  0.5- 1.0  Hz  and  1-2  Hz 
bands.  Symbols  are  the  same  as  in  Figure  3. 
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Figure  10.  Same  as  Figure  9  for  the  2-4  Hz  and  4-8  Hz  bands. 


Variations  in  Pg/Lg  with  Path 


One  of  the  inherent  difficulties  in  a  discrimination  study  is  accounting  for  path  effects. 
Since  the  earthquakes  and  explosions  are  very  often  located  in  different  regions,  it  is  even 
conceivable  that  observed  differences  in  phase  ratios  between  the  two  populations  may  be  due 
largely  to  differences  in  path  rather  than  differences  in  source  type.  The  logarithms  of  the 
uncensored  phase  ratio  values  for  four  frequency  bands  are  shown  on  a  map  of  the 
Califomia-Nevada  region  in  Figures  11  and  12.  The  mean  log  ratio  has  been  removed  to 
emphasize  variations.  Due  to  the  close  spacing  at  NTS,  the  mean  log  ratio  for  the  explosions 
is  shown  rather  than  the  individual  values. 

Two  main  groups  can  be  defined  according  to  the  source-station  travel  path:  a  northerly 
group  (above  36.5°  N),  which  includes  the  explosions  at  the  Nevada  Test  Site;  and  a  more 
southerly  group  consisting  of  earthquakes  along  the  Sierra-Nevada  fault  and  in  the  Mojave 
Desert.  In  the  lowest  frequency  band,  events  in  the  southerly  group  have  mostly  lower  than 
average  Pg/Lg  ratios.  In  the  1-2  Hz  band,  a  trend  of  increasing  Pg/Lg  with  distance  is  visible 
for  the  southerly  group,  which  grows  stronger  with  increasing  frequency.  Thus,  Pg  pro¬ 
pagates  more  efficiently  than  Lg  for  this  group  of  events.  This  discrepancy  may  be  due  either 
to  differences  in  geometric  spreading  (elastic  effects)  or  attenuation  (anelastic  effects).  The 
increase  of  the  discrepancy  with  increasing  frequency  suggests  that  attenuation  is  the  main 
cause.  On  the  other  hand,  sources  in  the  the  northerly  group,  including  the  explosions  at 
NTS,  have  almost  uniformly  higher  than  average  Pg/Lg  ratios  at  all  frequencies,  with  little 
apparent  distance  trend. 

The  most  striking  feature  of  these  maps  is  the  differences  between  the  two  groups.  The 
implication  is  that  the  path  exerts  a  dominant  effect  on  the  Pg/Lg  ratio.  Despite  the  fact  that 
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the  two  groups  of  sources  are  not  very  far  apart,  the  paths  from  source  to  receiver  are  quite 
different  due  to  the  Sierra-Nevada  batholith.  It  would  seem  that  the  Lg  wavetrain  for  the 
southern  group  should  travel  along  either  the  Sierra  Nevada  batholith  or,  more  likely,  along 
the  Great  Valley  for  most  of  its  path.  The  Lg  wavetrain  from  sources  in  the  northern  group, 
however,  first  traverses  basin  and  range  terrain,  and  then  the  Sierra  Nevada  batholith.  In 
essence,  Lg  from  the  first  group  travels  with  the  tectonic  and  topographic  grain,  while  Lg  from 
the  northern  sources  travels  across  the  grain.  Since  Lg  is  essentially  a  guided  wave,  the  latter 
path  type  is  likely  to  disrupt  Lg,  which  was  demonstrated  by  Kennett  et  al.  (1990)  using  3-D 
raytracing.  In  addition,  Lg  from  many  of  the  northern  sources  travels  through  the  vicinity  of 
the  Mono  Lakes  region,  where  high  attenuation  may  further  decrease  the  Lg  amplitude  relative 
to  Pg. 

In  contrast  to  the  differences  between  the  northern  and  southern  groups  of  events,  the 
average  explosion  Pg/Lg  is  actually  quite  similar  to  that  of  the  earthquakes  in  the  same  region 
in  most  of  the  frequency  bands.  Only  in  the  highest  frequency  band  does  the  average  Pg/Lg 
ratio  appear  to  be  higher  than  most  of  the  earthquake  Pg/Lg  ratios  in  the  same  region. 

We  have  computed  maximum-likelihood  averages  for  the  explosions  and  both  the  north¬ 
ern  and  southern  groups  of  earthquakes  (Figure  13,  Table  3),  using  the  technique  outlined  by 
Jih  and  Shumway  (1989).  The  error  bars  indicate  the  standard  deviations  of  the  maximum- 
likelihood  estimates.  Up  to  4  Hz,  the  explosion  population  overlaps  significantly  with  the 
nonhem  group  of  earthquakes,  while  it  appears  to  be  distinct  from  the  southern  group.  Note 
also  that  the  differences  between  the  northern  and  southern  groups  of  earthquakes  are  compar¬ 
able  to  those  between  the  explosions  and  northern  earthquakes.  This  suggests  that  the  travel 
path  is  as  important  as  the  source  type  in  determining  the  Pg/Lg  ratio.  In  the  highest  fre- 
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quency  band,  however,  the  two  earthquake  groups  converge,  and  the  explosion  Pg/Lg  ratio 
diverges  from  that  of  the  earthquakes. 


Figure  11.  Variations  in  log (Pg/Lg)  with  location  for  the  0.5-1.0  Hz  and  1-2  Hz  bands.  The 
mean  has  been  removed,  so  that  diamonds  indicate  lower  than  average  Pg/Lg  ratios  and 
crosses  indicate  higher  than  average  Pg/Lg  ratios.  For  clarity,  the  mean  of  the  explosions  is 
shown,  labeled  NTS,  instead  of  the  individual  values. 
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Figure  12.  Same  as  Figure  1 1  for  the  2-4  Hz  and  4-8  Hz  bands. 
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Frequency  Band 

Figure  13.  Maximum-likelihood  estimation  (MLE)  of  average  values  of  log (Pg/Lg)  for  three 
populations:  explosions  (asterisks),  earthquakes  from  the  southern  (California)  group  (dia¬ 
monds),  and  earthquakes  from  the  northern  (Nevada)  group  (triangles).  The  error  bars  indi¬ 
cate  the  MLE  standard  deviations. 
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Comparison  of  Explosions  and  Collapses 

One  case  does  exist  that  can  be  used  to  remove  the  path  effect  entirely:  a  comparison  of 
explosions  with  their  resultant  collapses.  Indeed,  the  only  difference  between  the  explosion 
and  the  collapse  should  be  the  source  function.  An  explosion  is  roughly  isotropic,  with  a  very 
short  time  history.  The  source  mechanism  of  a  collapse,  however,  is  less  well  known,  and 
may  range  anywhere  from  a  downward  point  force  at  the  cavity  to  a  spall-like  failure  of 
near-surface  rock  layers.  Thus  the  effective  depth  of  the  collapse  will  be  less  than  or  equal  to 
the  depth  of  the  explosion.  The  symmetry  of  the  collapse  source  is  most  probably  cylindrical 
rather  than  spherical  and  may  thus  be  expected  to  generate  significant  SV  waves. 

Three  of  the  explosions  in  this  data  set  had  collapses  which  were  well-recorded  at  JAS. 
A  comparison  of  the  Pg  spectra  for  the  TAJO  explosion  and  resultant  collapse  is  shown  in 
Figure  14.  The  spectra  have  been  smoothed  with  a  three-point  running  average.  The  collapse 
is  significantly  depleted  in  high  frequencies  relative  to  the  explosion,  to  the  point  where  the 
signal-to-noise  ratio  is  less  than  1  above  4  Hz.  Figure  15  shows  the  Pg/Lg  ratios  for  the  three 
explosion/collapse  pairs,  which  have  been  smoothed  with  a  five-point  running  average  and 
offset  by  a  factor  of  10  for  each  pair. 

The  collapse  ratios  track  those  of  the  explosions  fairly  closely  up  to  4  Hz,  despite 
presumably  disparate  source  functions,  emphasizing  the  dominant  roles  of  path  and  perhaps 
source  depth  in  determining  Pg/Lg  ratios.  Interestingly,  all  of  the  collapse  Pg/Lg  ratios  are 
lower  than  the  corresponding  explosion  ratios  at  the  low-frequency  end.  Recalling  that  explo¬ 
sion  Pg/Lg  generally  increases  with  depth  (Figures  7  and  8),  this  may  reflect  a  shallower 
effective  depth  for  the  low-frequency  component  of  the  collapse,  with  the  high-frequency 
component  having  a  depth  similar  to  that  of  the  explosion. 
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Figure  14.  Comparison  of  TAJO  collapse  (heavy  line)  and  explosion  (thin  line)  Pg  spectra. 
The  dashed  line  is  the  noise  spectrum  of  the  collapse  seismogram.  All  of  the  spectra  have 
been  smoothed  with  a  three-point  running  average. 
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Figure  15.  Comparison  of  explosion  (solid  line)  and  collapse  (dashed  line)  Pg/Lg  ratios  as  a 
function  of  frequency.  BRETON  and  TAJO  have  been  offset  by  respective  factors  of  10  and 
100  from  DUORO.  All  ratios  have  been  smoothed  with  a  five-point  running  average. 


CONCLUSION 


The  results  above  shed  some  light  on  why  the  performance  of  the  Pg/Lg  disci iminant  is 
so  uneven.  In  this  case,  sources  were  selected  to  be  as  close  to  each  other  as  possible,  thus 
reducing  factors  induced  by  azimuthally  varying  receiver  functions  and  ideally  mitigating  the 
effects  of  the  travel  path.  Nevertheless,  sources  located  as  little  as  100  km  apart  have 
significantly  different  Pg/Lg  ratios.  If  the  set  of  sources  is  restricted  only  to  sources  that 
travel  roughly  the  same  path  (the  northern  group  of  earthquakes),  then  the  differences  between 
earthquakes  and  explosions  decrease  dramatically. 

These  results  also  present  grounds  for  both  skepticism  and  hope  in  Pg/Lg  discriminants. 
On  one  hand,  the  Pg/Lg  ratios  for  earthquakes  and  explosions  seem  to  be  controlled  more 
strongly  by  the  travel  path  than  the  source  type  at  most  frequencies.  Complicating  matters 
further  is  the  dependence  of  the  explosion  Pg/Lg  ratio  on  source  depth  and/or  magnitude, 
which  tends  to  make  shallower,  smaller  explosions  look  more  like  earthquakes.  On  the  other 
hand,  the  highest  frequency  range  studied,  4-8  Hz,  seems  to  show  some  promise  in  that  the 
Pg/Lg  ratio  depends  somewhat  less  on  path  than  in  other  frequency  bands.  Also,  the  explo¬ 
sion  population  is  better  separated  from  the  earthquakes  in  this  band.  Still  higher  frequencies 
may  improve  this  separation,  but  they  could  not  be  used  for  JAS,  for  which  the  sample  rate 
was  only  20  samples/second.  However,  the  Pg/Lg  ratio  at  higher  frequencies  also  shows  a 
stronger  dependency  on  distance,  probably  due  to  the  differential  attenuation  of  the  two 
phases.  Thus,  distance  corrections  may  further  increase  the  separation  of  earthquakes  and 
explosions,  with  the  caveat  that  these  distance  corrections  will  likely  be  path-dependent. 
Another  drawback  to  the  higher  frequencies  is  an  apparent  increase  in  scatter  for  both  earth¬ 
quakes  and  explosions,  and  in  some  cases  a  lower  signal-to-noise  ratio. 
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